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We examine the inference of quantum density operators from incomplete information by means 
of the maximization of general non-additive entropic forms. Extended thermodynamic relations are 
given. When applied to a bipartite spin | system, the formalism allows to avoid fake entanglement 
for data based on the Bell-CHSH observable, and, in general, on any set of Bell constraints. 
Particular results obtained with the Tsallis entropy and with an introduced exponential entropic 
form are also discussed. 
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The relation between two fundamental concepts, en¬ 
tropy and quantum entanglement, has recently aroused 
great interest in quantum information theory [H-Q. A 
system composed of two subsystems A and B is said to 
be unentangled or separable , if and only if the density op¬ 
erator p can be written as a convex combination of uncor¬ 
related densities, i.e., p = qjpf ® pf , with qj > 0. In 
this case the system admits a local description in terms 
of hidden variables. Otherwise, it is said to be entan¬ 
gled or inseparable. The system becomes then suitable, 
in principle, for applications like quantum cryptography 
[9j and teleportation [10]. 

When the available information about the system is in¬ 
complete, consisting for instance of the expectation values 
of a reduced set of observables, one faces the problem of 
first determining if entanglement is actually implied by 
the data, and then selecting the most probable or rep¬ 
resentative density operator compatible with these data. 
An ideal inference scheme in this scenario should then 
a) avoid fake entanglement [lj, i.e., should not yield an 
entangled density if there is a separable density that re¬ 
produces the data, b) be least biased , in the sense that 
some measure of lack of information is maximized, and 
c) be simple enough to be readily applied. As shown 
in (lj, the standard approach based on the direct max¬ 
imization of the von Neumann entropy S = —Trplnp, 
does not comply with a) already for two spin i systems. 
The essential reason is that this entropy is not a good en¬ 
tanglement indicator 0, §j, even in those cases where 
entanglement is fully determined by the eigenvalues of p. 
A solution was also provided in [lj: one should first de¬ 
termine the set of densities that minimize entanglement, 
and then maximize entropy within this set. Although 
certainly rigorous, this procedure is not easy to imple¬ 
ment in general, and departs conceptually from a more 
basic approach based on the maximization of a single 
information measure. 

As is well known, the von Neumann entropy is based 
on the Shannon information measure, which is the unique 
one satisfying the four Khinchin axioms ;n]- However, if 
the fourth axiom, which is concerned with additivity, is 
lifted, other information measures become feasible. The 
most famous recent example is the Tsallis entropy 0, 


which has been applied to a wide range of phenomena 
characterized by non-extensivity [ju|, including recently 
the problem of quantum entanglement Si- 

The aim of this work is first to discuss more general 
non-additive entropic forms, based on arbitrary concave 
functions, and the ensuing densities that maximize these 
forms subject to given constraints. Though sharing many 
features with the usual von Neumann based statistics, 
the extended formalism opens new possibilities, in par¬ 
ticular that of approaching, for certain functions, den¬ 
sities whose largest eigenvalue has the minimum value 
compatible with the available data. For a system of two 
qubits, this allows to satisfy previous requirements a), 
b), c) simultaneously for any set of Bell constraints, by 
means of a single maximization. In particular, for infor¬ 
mation based on the Bell-CHSH observable [l|, we will 
show that fake entanglement can be averted even with¬ 
out including data about the dispersion (in contrast with 
H HI), for a wide class of entropic functions. As par¬ 
ticular cases, we will examine results obtained with the 
Tsallis entropy, applied here with standard expectation 
values (as opposed to Q, where the g-expectation values 
were used), and will also introduce an exponential en¬ 
tropic function, which will provide fully analytic results. 

Given a density operator p = YliPi N)(*Ij with Vi > 0, 
Tr p = Y^iPi = 1 and the sum running over a complete 
set of orthonormal states QA |*)(z| = /), we will consider 
the entropic form 

Sf{p) = Tr /(p) = X /(p*) ’ W 

i 

where / is a smooth concave function < 0 for 

p € (0,1) and f(p) continuous for p £ [0,1]) satisfy¬ 
ing /(0) = /(1) = 0. Although Eq. (JT|) is certainly 
not the most general information measure, it is the 
most simple generalization of the von Neumann entropy 
(/(p) = — plnp) and comprises useful extensions. With 
the exception of additivity, the basic properties of en¬ 
tropy are satisfied for an arbitrary function / of this 
form. In particular, Sf(p) > 0, with Sf{p) = 0 only 
in the case of a pure state ( p 2 = p), its maximum is 
attained for a uniform distribution (pi = 1/ra, with n 
the space dimension), with the maximum value nf(l/n) 
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increasing with n , and it is not affected by states with 
Pi = 0. Concavity of / implies concavity of Sf(p) (l4| . 
i.e., S f Q2j QjPj) > Ej ). with Qj > E, % = 1, 

as well as the important property that Sf(p) cannot de¬ 
crease by dephasing , i.e., by removing the off-diagonal 
elements of p in any basis of orthonormal states | i a ): 

Sf(p) < S f (p 0 ), p 0 = J2(io\p\io)\io){io\- (2) 

i 

A sufficient condition for Sf to be sub (super) additive 
is that pf"(p) be a decreasing (increasing) function of p 
for p G (0,1), since in this case f(pq) — qf(p)—pf(q) < 0 
(> 0) V p, q G [0,1], implying [H 

Sf(.P A ®P B ) > Sf( P A ) +Sf(p B ) if (pf"(p))' >0. (3) 

Additivity amongst entropies of the form ( 1 ) holds only 
if (p/"(p))' = 0, which leads immediately to the Shannon 
form f(p) = —kplnp, k > 0 . 

Maximization of Sf(p) subject to the constraints of 
to + 1 expectation values (O a ) = Tr pO a , a = 0,..., to, 
where O a are linearly independent observables (not nec¬ 
essarily commuting) and we have set Oq = I to account 
for normalization ((/) = 1 ), leads to the maximization of 


respectively, we obtain the thermodynamic relationships 


dSf(X) 
dX a 
dSfjO ) 
d(O a ) 


= ~(O a ) , 


= A n 


d 2 S f ( A) 

dX a dX/s 


( 7 ) 


d 2 S f (0) 

d(o a )d(Op) 


8 ) 


AolP — \Oq, \jh) (jh \ 0p\ih)Cij ? (9) 

id 


Cij = -5 ijP \hi) + (1 - > 0 • ( 10 ) 


Only the second derivatives in © © depend explicitly 
on /, through the first term in (flOl) . As C l; j > 0, A aa > 
0. Moreover, all eigenvalues of A are non-negative , i.e., 
A(x y | {ih\Qa I jh .) | with Q a UBgO B and 

U defined by \U tr AU] a p = A a 5 a p. Hence Sf is a convex 
function of A, whereas Sf a concave function of O, as in 
the von Neumann case. If [O a ,Op\ = 0 V a,/3, Eq. © 
becomes A a p = -Tr p' O a O 0 , with p' = EiP'(hi)\ih)(ih\ 
(for f(p) = —p Inp, p'[hi) = -pi and p' = -p). 

We will be here interested in functions of the form 


f[p) = k[p-g q [p)), ( 11 ) 

where k > 0 and g q (p) is a convex [g q [p) > 0 ) increasing 
function satisfying g q [ 0 ) = 0 , g q [ 1 ) = 1 , and 


s f[p) = S f [p) - ^2 X a (O a ) = Tr [f[p) - ph], (4) Xm^g^Pi)/g^pj) = 0 if p t < p 3 . (12) 

a 

i 7 tx r -i. • i • Hence, for sufficiently large q (and finite dimension n) 

where ri = and A a are Lagrange multipliers. 

Writing h = hi\i h ){i h |, the maximum is attained for 5/(p) = jfe[l - Tr g q [p)} rs k{ 1 - n M g q [PM)) , 


p = p[h) = ^2pi\ih){ih\ , 
i 


Pi = p[hi) = 


\fV\hi) 

0 


/'(l) <hi< /'(0) 
hi > /'( 0 ) 


( 5 ) 

( 6 ) 


where [/ , ]” 1 is tf 16 inverse of the function /'. The cut¬ 
off for hi > /'( 0 ) obviously arises only when /'( 0 ) is 
finite , and is the main difference with the von Neumann 
case (where Eq. © becomes the exponential distribution 
Pi = e~( hi+1 \ with hi > —1). Nevertheless, due to the 
concavity of f,Pi is always a non-increasing function of 
K {p' (hi) = 1/ f"(Pi) < 0 if /'(1) < hi < /'(0), and 0 if 
hi > f'( 0)) that vanishes for hi —> oo V /. 

Eqs. © © can be easily derived. As Tr ph = Tr p^h, 
with p h = ESh\p\ih)\ih)(ih\, Eq- © implies S f (p) < 
Sf(ph)- The optimum density satisfies then [p, h] = 0, so 
that p h = p and Sj(p) = EifiPi ) - Pihi- Due to the 
concavity of /, this will have a unique maximum for pi G 
[ 0 , 1 ], determined by f(pi) = hi if /'( 1 ) < hi < f'( 0 ), 
or otherwise located at one of the borders, which leads 
to Eq. © (for a non-standard normalization (I) > 1, an 
upper cutoff pi = 1 if hi < f'(l) would also apply). 

Considering Sf(p) and Sf(p) at the maximum © as 
functions of = (A 0 ,...,A m ) and 6 = ((O 0 ),..., (O m )) 


where pm is the largest eigenvalue of p and % its mul¬ 
tiplicity. The density that maximizes Sf(p) (i.e., mini¬ 
mizes Tr g q (p)) subject to a given set of constraints, will 
then possess, for q —> oo, the minimum largest eigen¬ 
value Pm compatible with the available data, as in this 
case nMg q (PM) is minimum. This property may in fact 
be fulfilled already for finite values of q (i.e., typically 
q > q c ) in some cases, as will be seen below. 

Similarly, maximization of the entropy associated with 

Kp) =/(1-p) = fc(l-p-Sg(l-p)), (13) 

which is also concave and satisfies /(0) = /( 1 ) = 0, leads 
to a density which possesses, for q —> oo, the maximum 
smallest eigenvalue compatible with the available data. 
In this limit g q ( 1 — Pi)/g q (l — Pj) —> 0 if pi > pj, so that 
Sj(p) « k[n — 1 — n m g q ( 1 - p m )\ for large q, with p m 
the smallest eigenvalue and n m its multiplicity. This is 
maximum if p m is maximum. 

As a particular example, we have in first term 

f(p) = (p-p 9 )/(q- 1), 9 > 0 , (14) 

which is concave for q > 0 , approaches —plnp for q —> 1 
and is of the form (HID for q > 1, satisfying (fT21) . In 
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this case, Sf(p) = (1 — Tr p q )/(q — 1) is the Tsallis en¬ 
tropy, which is sub (super)-additive for q > 1 (q < 1 ), 
in agreement with Eq. © (( pf"{p)Y = q (i - q)p q 2 )- 
The q = 2 case is particularly relevant, since maximiza¬ 
tion of Sf(p) becomes equivalent to the minimization of 
Trp 2 = (^T <.(pi — pj) 2 + 1 )/n, i.e., to a least squares 
condition for the probabilities or their differences. For 
q = 2, Sf(p) also coincides with the information measure 
of ref. fg|. Eq. © becomes, setting h c = f'{ 0), 

Pi = {[1 ~{q~ 1 )hi]/q} 1/i - q ~ 1) , -1 <hi<h c 


which determines p+, and hence p- = — b, po = ^(1 + 

b) — p+, for a given /. If /'(0) is finite, a root of Eq. (fl9l) 
for e [ 6 , A (1 + 6 )] will exist only if | 6 | < b c , with b c 
the root of 

/ , (M + /'(0)-2/'(^^)=0. (20) 

Eq. (l20l) is equivalent to /'(0) = Ao — Ai, and determines 
the onset of the cut-off for . It has a single root b c £ 
[|, 1], with 6 C —>• 1 if /'(0) —> oo. For b > b c , A 0 — Ai > 
/'( 0 ) and we obtain the solution 


and pi = 0 if hi > h c , with h c = (oo) if q > 1 (q < 1 ). 

Another example is the exponential function 

f{p) = q~ l [p- (e 9P - l)/(e 9 - i)], (15) 

which is concave for any q , satisfies /( 0 ) = /( 1 ) = 0 and 
approaches \p{ 1 — p) for q —> 0, i.e., proportional to the 
Tsallis case q = 2. For q > 0, it is of the form HID 
and fulfills Eq. (fl2l) . Moreover, f~ q (p) = f q ( 1 — p). As 
( pf"(p))' = e qp (l -1- qp)q/(l — e q ), Sf(p) is sub-additive 
for q > — 1. Eq. © becomes 

Pi = q~ x ln[l - {e q - 1 )(/i, - 6 C )], -1 < hi - h c < 0 

and Pi = 0 if hi > h c , with h c = | > 0 . 

Application. Let us consider now a bipartite spin A 
system. Starting from the basic unentangled states |tt), 
I tl) > | It) 5 | It) , the Bell basis is formed by the fully 
entangled orthonormal states |\E f=fc ) = (|tl) ± \i'\))/\/2, 
|$ ± ) = (|tt) ± |||))/\/2 (|T-) is the singlet while |4-+), 
l*^) are spin 1 states with (S) = 0). We will consider 
here Bell constraints [l|, i.e., mean values of observables 
which are diagonal in the Bell basis. Let us first examine 
the case of ref. flj , where the available information is the 
expectation value of the (scaled) Bell-CHSH observable 

B=|$ + )<$+|-|W-)<^-|. (16) 

According to Eq. ©, the density p that satisfies 


Trp = 1, Tr pB = b, | 6 | < 1, (17) 


and maximizes © is of the form 

p = p{\oI + A rf) = po(|'F + )('F + | + |$-)($-|) + 

p + |$+)($+|+p_|^-)(^-|,(18) 

where po = p{ Ao), p± = p{ Ao ± Ai). The constraints (1171) 
become just 2po+p+ +p~ = 1, p+—p- = b. We may 
consider b > 0, in which case Ai < 0 and p+ > po > P-, 
since for b —> — b , Ai — > —A± and p± —> p=p. 

If f'( 1) < Ao ± Ai < /'(0), there is no cut-off and the 
constraints lead to the single equation 


f\p + ) + f\p + -b)-2f\ 


1 + b 


P+) = 0, |&| < b c (19) 


p+=b, p-= 0 , p 0 = \(l-b), b c <b< 1 . ( 21 ) 

Eqs. (flTH) (T5T1) become apparent from the entropy 

Sf(p) = f(?+) + f{P+ ~b)+ 2 /(-y- - P+) ■ (22) 

For fixed b > 0, Eq. (B^l) is a concave function of p + for 
p+ £ [b, ^(1 + &)], with its maximum located within the 
interval if | 6 | < b c , being then determined by (flfll) . and 
at the left border if b c < b < 1 , leading to m- At the 
maximum, Ai = dSf(p)/db = f'(p +) — f'(po ) in both 
cases, with Ao = f'(po)- Eqs. (fTUl) (Ell) imply that p + is 
an increasing function of b for b > 0 . 

For b -£ 0, Eq. (HU) leads to 

p + = A(l + 26 + 7 5 2 ) + 0 ( 6 4 ), 7 => ( 23 ) 

where 7 < 1 (> 1 ) if Sf is sub (super) additive and 
satisfies Eq. ©■ Hence, for b = 0 we obtain the uniform 
distribution p + = = po = 7 for any /. On the other 

hand, for b —> 1 , p + —> 1 and p -» |$ + )($ + |. 

The important question that now arises is whether for 
a given /, the previous scheme gives fake entanglement. 
A density p diagonal in the Bell basis is unentangled if 
and only if its largest eigenvalue is less than ~ Q. The 
density of the form (1181) that complies with (tT71) and pos¬ 
sesses the minimum largest eigenvalue corresponds to 

P+ = Po = |(1 + b ), p _ = 1(1 - 36), 0 < b < A . 

P+ = b, p- = 0, po = A(i - b ), g < b < V 

where > po > p~. Unentangled solutions are then 
feasible only if b < A. Note also that for Bell constraints, 
entanglement is minimized by densities which are diago¬ 
nal in the Bell basis [lj, so that no unentangled density of 
any form complying with C3 exists for 6 > A. It is now 
seen from m that when /'( 0 ) is finite, the maximum en¬ 
tropy density coincides with (El for b > b c > A. Hence, 
as is an increasing function of b , fake entanglement 
will be avoided for those f for which b c < ^. 

Particular solutions. In the von Neumann case, Eq. 
(fl9l) yields p + = y(l + b) 2 , with b c = 1, in agreement 
with ref. lj and Eq. (E^l) (7 = 1). Fake entanglement 
occurs for \/2 — 1 < b < \. 
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In the Tsallis case (fl4l) . /'(0) is finite for q > 1 and 
Eq. (l20l) leads to 

b c = [ 1 + 2 1 - 1 /^- 1 )]- 1 , q> 1, (25) 

which is a decreasing function of q satisfying b c < A for 
q > 2. Hence, fake entanglement will be avoided V q > 2. 
For q = 2, the solution of Eq. m is specially simple, 

P+ = 1(1 + 26), 0 < b < i, p+=b, \<b< 1 (26) 

which is in agreement with (1231) (7 = 2 — q) and rep¬ 
resents the solution of minimum squares. The onset of 
entanglement occurs here exactly at b = b c . 

Although a simple analytic solution of (flT)l) for arbi¬ 
trary q is not feasible, it is easy to verify that Eq. (PI is 
obtained for q —> 00 V b. In this limit, b c —> while Eq. 

(flfll) yields, for large q, p+ ~ |(1 + 6)[1 + 
which approaches |(1 + b) for q — > 00 . 

For the exponential function (fT5l) . the solution of Eq. 
dm is analytic for any q , 

P+ = |(1 + 2b) — — In cosh (^bq) , 0 < b < b c (27) 
Zq 

b c = \ + - q HP q - e - 3^-]. P q = [i + \ji + ^] 1/3 , 

with p + = b for b > b c . For q —> 00, b c —> A and J27} leads 
immediately to the solution with the minimum largest 
eigenvalue, Eq. (1241) . Again, b c is a decreasing function 
of q , with fe c < A for q > 0, so that fake entanglement 
is here avoided V q > 0. For q —> 0, b c —> A and Eq. 
(1271) reduces to G51) . For b —> 0, p + ~ |(1 + 2b — \qb 2 ), 
in agreement with (1251) . Finally, for q —> —00, b c —»• 1 
and p+ —»• 4(1 + 36), with p_ = po = A(1 — 6), which is 
the solution with the maximum smallest eigenvalue of p. 
This gives fake entanglement for b £ [A, ^], the maximum 
interval for maximum entropy densities, as p+ is in this 
case maximum. 

Inclusion of the dispersion. If the dispersion of B 
is also provided Q, through the expectation value of 
B 2 = |<f >+ )($ + | + the final maximum entropy 

density is actually independent of the choice of f. In this 
case p = p( Ao + Ai B + X 2 B 2 ) is also of the form (1T51) . 
with p± = p( Ao ± Ai + A2), Po = p{ Ao), which are com¬ 
pletely determined by the constraints, i.e., p± = A(6 2 ±6), 
Po = 5(1—62), where b 2 = Tr pB 2 = p + +p_. The only role 
played here by maximum entropy is to impose a density 
of the form (II 8 [I . which holds for any /, and fake entan¬ 
glement is then always avoided. Note also that when only 
b is given, the solution m implies minimum dispersion , 
as b 2 = 2p- + b is minimum (see also ref. |). 

General Bell constraints. Higher values of q , or even 
the limit q —> 00, may be required in general to avoid 
fake entanglement. For example, if B a = |<I> + )(<1> + | — 
ex 141 ) (\k |, with a > 0 and Tr pB a = b> 0, unentangled 
densities complying with these data are again feasible 


only if b < ^. However, in the maximum entropy density, 
b c will increase as a decreases below 1, with b c —> 1 
for a —> 0. This implies a higher threshold to avoid 
fake entanglement, namely q > 1 + log 2 (l + a -1 ) in the 
Tsallis case (THl) and q > —41n(a) in (fl5l) . For a —>■ 
0, both values become 00, although the interval of fake 
entanglement vanishes, approaching [A — Aa, if b c > A 
(for a = 0, p+ = b, p- = po = ^(1 — b) V /, becoming p 
a Werner state, and fake entanglement does not occur). 


The present arguments are valid for any type of Bell 
constraints. In this case, densities of minimum entan- 
glement, as measured by the entanglement of formation 
Ef{p) [ijl, are diagonal in the Bell basis l|, and possess 
the minimum largest eigenvalue pm compatible with the 
available data if pm > y ■ This is so because Ep{p ) is an 
increasing function of the concurrence C(p), which for a 
system of two qubits reduces to 2 pm — 1 if pm > 4 (and 
0 otherwise) when p is diagonal in the Bell basis [5|, 17] - 
Maximum entropy densities constructed with functions 
satisfying G3J HU]) will then possess minimum entangle¬ 
ment for q —> 00 , although in some cases this may hold 
already for finite values of q , as seen in the example. For 
sufficiently large q, the entropies Sf(p) will be essentially 
decreasing functions of pm, being then good entanglement 
indicators for these densities. This may also be the case 
in systems of n qubits for special constraints that lead to 
densities diagonal in a basis of fully entangled states (like 
the GHZ states used in Iff), where separability is again 
favored by low values of the largest eigenvalue. Further 
investigations in this direction are in progress. 


In summary, we have examined the use of general non¬ 
additive entropic forms for the inference of quantum den¬ 
sity operators. The formalism enables a direct way to 
infer least biased densities with minimum entanglement 
for a system of two qubits, and hence to avert fake en¬ 
tanglement, when the information consists of any set of 
Bell-constraints. It also provides a general framework for 
the description of the thermodynamic aspects of entan¬ 
glement, as well as a more deep foundation of the success 
that non-additive entropies like that of Tsallis may en¬ 
counter in this type of problems. 
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